Optical effects of spin currents in semiconductors 
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A spin current has novel linear and second-order nonlinear optical effects due to its symmetry properties. 
With the symmetry analysis and the eight-band microscopic calculation we have systematically investigated 
the interaction between a spin current and a polarized light beam (or the "photon spin current") in direct-gap 
semiconductors. This interaction is rooted in the intrinsic spin-orbit coupling in valence bands and does not 
rely on the Rashba or Dresselhaus effect. The light-spin current interaction results in an optical birefringence 
effect of the spin current. The symmetry analysis indicates that in a semiconductor with inversion symmetry, the 
linear birefringence effect vanishes and only the circular birefringence effect exists. The circular birefringence 
effect is similar to the Faraday rotation in magneto-optics but involves no net magnetization nor breaking the 
time-reversal symmetry. Moreover, a spin current can induce the second-order nonlinear optical processes due 
to the inversion- symmetry breaking. These findings form a basis of measuring a pure spin current where and 
when it flows with the standard optical spectroscopy, which may provide a toolbox to explore a wealth of physics 
connecting the spintronics and photonics. 

PACS numbers: 72.25.Dc, 78.20.Ls, 42.65 An 78.30.Fs 



I. INTRODUCTION 

A pure spin current consists of flows of opposite spins in 
opposite directions with the same amplitude. It bears neither 
net charge current nor net spin polarization. Spin currents are 
a key element in spintronics.^ ^ 

Detection of spin currents is important for characteriza- 
tion and applications in future spintronics technologies.^^ 
While a polarized spin current may be detected by the con- 
ventional Faraday/Kerr rotation spectroscopy^ ^ or through 
ferromagnetic filters!^ a pure spin current, without a di- 
rect electromagnetic induction, is much less traceable. Still, 
pure spin currents have been detected in a few pioneering 
experiments in which they were converted into signals de- 
tectable by conventional techniques. For example, the spin- 
polarized electrons or excitons accumulated at the sample 
edges where a spin current is terminated may be detected 
by the Faraday/Kerr rotation,'^ polarized light emission,^ 
and polarization-selective absorption. Or the inverse spin 
Hall effec i^-^tti6| ^^^^ covert a spin current into 

charge/voltage signals for electric measurement.'^^''^ All such 
measurements, however, disturb the spin currents to some 
extent and are indirect. We are motivated to find a non- 
destructive way to directly measure a pure spin current.^^^^ 

A basic symmetry principle states that whenever there is a 
current breaking the fundamental symmetries of a system, an 
interaction may arise between the current and another current 
of the same symmetry-breaking type so that the fundamental 
symmetries are retained.^^ A classic example is the Ampere 
effect and the 0rsted effect where a charge current is coupled 
to another charge current or a magnet. A straightforward ana- 
logue suggests that a pure spin current may be coupled to 
another spin current. Such an idea stimulated the proposal 
of direct measurement of a pure spin current in a direct-gap 
semiconductor by a polarized light beam.^^ A polarized light 



beam can be regarded as such a "photon spin current'!^ by 
mapping the photon polarization into a spin- 1/2 in the Jones 
vector representatioil^ 

cos ^e'^^^^n^ + sin ^e-'^^^^n. ~ |6>, 0) , (1) 

where the right/left circular polarization n+/_ corresponds to 
the spin up/down state IT / vl) quantized along the light propa- 
gation direction. The effective interaction between a pure spin 
current and a polarized light causes a phase delay which de- 
pends on the light polarization and wavevector. The observ- 
able result is a circular birefringence effect which is similar 
to the Faraday rotation but involves no net magnetization nor 
time-reversal symmetry breaking. Since Faraday's disocvery 
in 1845,^^ the circular birefringence effect of spin currents 
is the first example of Faraday rotation without time-reversal 
symmetry breaking. We dub this effect as spin current Fara- 
day effect. 

At this point, we should mention a recent remarkable ex- 
periment realizing the direct in-situ detection of a spin current 
through the Doppler effect of a spin- wave .'^ In fact, the ob- 
served Doppler effect and our predicted optical birefringence 
effect are fundamentally related to each other. The former is 
the frequency shift of the spin wave, while the latter is the 
phase shift of the light accumulated by a frequency shift over 
a coupling time. The frequency shift is measured in the near- 
field as in the experiment, while the phase shift should be 
measured in the far field by light polarization detection. Fun- 
damentally, both are due to the effective coupling between a 
pure spin current and another "probe spin current", either a 
spin wave or a polarized light, mediated by virtual excitations 
in the systems. 

The effective light- spin current interaction is induced in a 
semiconductor by virtual excitations of electron-hole pairs. 
The specific form of the phenomenological coupling depends 
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on the microscopic mechanisms!^ Since the Hght polarization 
essentially couples only to the orbital motion of electrons, the 
spin-orbit interaction is needed to establish the effective cou- 
pling. As there is inherent spin-orbit coupling in the valence 
bands due to the relativity effect, the Rashba or Dresselhaus 
effect due to the spatial inversion asymmetrie^^^^ is not a 
necessity, thus the system can bear the inversion symmetry. 

The optical birefringence effect of spin currents^^ is usually 
very weak, because a tiny light wave vector q is involved in the 
coupling to the velocity v of the pure spin currents. However, 
if the velocity of spin currents couples to another optical field, 

q . V ^ F2 • V 

the coupling will be much enhanced. This means we can use 
the second optical field to drive the spins, which may result 
in the nonlinear optics of the pure spin current. In fact, such 
an analogy stimulated the prediction of the second-order non- 
linear optical effects of pure spin currents,^^ which was soon 
verified by experiments.^^ 

In Refs. 20 and 21 we have sketched the main ideas based 
on symmetry arguments and given the key expressions in a 
special model neglecting the energy band anisotropy. In this 
paper, we will investigate in a more comprehensive way the 
linear and second-order nonlinear optical effects of pure spin 
currents, including a systematic symmetry analysis of all rele- 
vant physical quantities, and a detailed derivation for the effec- 
tive Hamiltonian as well as the second-order nonlinear optical 
susceptibility. The microscopic derivation confirms the quali- 
tative results obtained by the symmetry analysis. In particular, 
both the symmetry analysis and the microscopic calculation 
lead to the conclusion that the linear birefringence effect (sim- 
ilar to the Voigt effect in magneto-optics) always vanishes and 
only the circular birefringence effect exists, and the energy 
band anisotropy induces only a relatively small quantitative 
modification of the results. The absence of the Voigt effect is 
fundamentally related to the lack of the |0) state in the physical 
spin of photons [not the pseudo-spin in Eq. ([T])]. The micro- 
scopic mechanism of both linear and second-order nonlinear 
effects can be understood in a unified physical picture. 

In this paper, we assume that the host semiconductor sys- 
tem has the inversion symmetry. We note that in compound 
semiconductors such as GaAs the inversion symmetry is bro- 
ken, which, though a small effect, is critical to some schemes 
of spin current generation. In our present scheme, how- 
ever, the small inversion asymmetry in compound semicon- 
ductors is not important. For conditions used in our micro- 
scopic calculation, the spin splitting resulting from the Dres- 
selhaus effect due to the bulk inversion asymmetry (The Dres- 
selhaus splitting is ~ 0.01 meV in GaAs with doping density 
~ 10^^ cm"^), much less than the detuning of the light from 
the interband transitions^^ that mediate the effective interac- 
tion, so we can neglect the bulk inversion asymmetry in the 
measurement process even though it could be of vital impor- 
tance in generating the spin current. Also, in this paper, we 
consider only bulk materials, so the structure inversion asym- 
metry plays no role, though it is the basis of the Rashba effect. 
Without considering the Dresselhaus and Rashba effects due 
to inversion asymmetries, we avoid the subtlety in the defini- 



tion of a spin currentl^^'^^'The effect of inversion asymmetries 
on the interaction between the polarized light beams and a 
spin current, of course, is worth further study, but we prefer 
leaving this question open in this paper. 

The paper is organized as follows. Sec. |ll| presents a sys- 
tematic symmetry analysis for the coupling system to give a 
qualitative understanding of the linear and circular birefrin- 
gence effects and the second-order nonlinear optical effect of 



pure spin currents. Sec. Ill gives the theoretical model and mi- 
croscopic derivations for both the linear and the second-order 
nonlinear optical effects, and also explains the physical pic- 
tures for the microscopic mechanism of optical effects of spin 
currents. Sec. [IV] presents the numerical results and discus- 



sions of the experiment scheme. Sec. [V| concludes this paper. 



II. SYMMETRY ANALYSIS 

We will particularly consider the time-reversal (T) and the 
space-inversion (^) symmetries of all the relevant physical 
quantities, and the geometry symmetry for a specific form of 
spin currents. According to the symmetry analysis, a pure 
spin current may result in a circular birefringence effect but 
not a linear birefringence effect, and as it breaks V symmetry, 
a spin current can induce the second-order nonlinear optical 
processes. 



A. Linear optical effects 

We assume the whole system has the T and P symmetries 
at equilibrium. Namely, the effective coupling between a spin 
polarization or a spin current in the semiconductor system and 
a probe should have both symmetries, i.e., the transformation 
properties of the effective Hamiltonian *Keff are 





T 


V 









(2) 



where -1-/- refers to even/odd under the corresponding sym- 
metry transformations. 

In our study, a pure spin current is made of a non- 
equilibrium distribution of spin polarization in the momentum 
space. In general, it can be quantified by a rank-2 pseudo- 
tensor defined by (with volume of the material taken as unity) 



(3) 



where Sp is the spin polarization and Vp is the velocity of a 
particle with wave vector p, and e is the electron charge. The 
"photon spin current" tensor for a polarized light beam with 
electric field F(r, t) = (F+n+ -h F_n_) + c.c. is formu- 

lated as 



I = Iq = ^ (^xXZ -h /yyz -h 1{LT^ , 



(4a) 
(4b) 
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where q is the wave vector of the Hght beam, the unit vector 
z is chosen along the direction of q so that q = qz, the unit 
vectors x and y are related to the light polarization through 
n+ = (=Fx - iy) I ^Jl, and cr^ (j = x, y, z) is the Pauli matrix. 
For completeness, we also consider the spin polarization of 
the system 

S = ^Sp. (5) 
p 

The transformation properties under the T and P of the rel- 
evant physical quantities are 
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It is worth mentioning here there is no |0) state in the physical 
spin of photons, and the photon pseudo spin and Iy do not 
break the 7~- symmetry, for it involves the 2nd-order spin flip 
processes such as | + 1) |0) I - 1). In the following 
we will use these quantities to form an efl'ective Hamiltonian 
(undetermined up to a few coupling constants) satisfying the 
T and P symmetries. Since the interaction of the light with a 
spin is usually weak, we only consider the efl'ect in the leading 
order, which is bilinear in the spin and light quantities. 

Net spin polarization. The only optical quantity of the same 
symmetry -breaking type as the spin polarization is I^z. Thus 
the efl'ective interaction between a spin polarization and a light 
beam has the form 

= fo/zS . z, (7) 

with a coupling constant to be determined by the specific 
microscopic mechanism. Such a coupling corresponds to the 
conventional Faraday eff'ect in magnetooptics.^^ We would 
like to point out here that a spin polarization could induce the 
Voigt eff'ect. In order to have the same symmetry-breaking 
type for and Iy, the spin polarization should be of an even 
power. Thus to the leading order, the eff'ective interaction has 
the form 

This explains that the Voigt eff'ect is quadratic in the spin po- 
larization or the applied external magnetic field. 

Pure spin current. There is no term in the light polarization 
Ij (j = X, y, z) that has the same symmetry -breaking type as 
the spin current, so it is not possible to have linear (nonlinear 
optics is of course possible) interaction between the spin cur- 
rent and the light without involving the wave vector. Consid- 
ering the wave vector of the light, coupling between the spin 
current and the photon current ql^zz is possible. The linear 
birefringence eff'ect (similar to the Voigt eff'ect in magnetoop- 
tics) is absent. Due to the lack of |0) state in the physical spin 
of photons, Ix and Iy preserve T symmetry. Therefore there is 
no linear coupling of J to ^/^xx and qlyjy. 

Furthermore, if the system has spherical symmetry, the ef- 
fective Hamiltonian would have a simple tensor contraction 
form as 

Hff = f i^^zTr (J) + • J • z, (9) 



with only two coupling constants f i and ^2 to be determined 
by the microscopic mechanisms. A possible spherically sym- 
metric system is the vacuum, but in general a semiconductor 
as a crystal does not have this symmetry. The general eff'ective 
interaction in a semiconductor should have the form 

"Hl^ = ^/,zz : ^ : J, (10) 

where is a parameter tensor determined by the microscopic 
structure of the material. Since only the light polarization term 
Iz appears in the interaction, the optical birefringence eff'ect is 
circular, similar to the Faraday rotation. 

In realistic case, the spin current often has some special 
form. As a general case a spin current tensor can have the 
form as 

J = JxXZ -h JyYZ -h JzZZ = JZ, (11) 

where Z is the unit vector along the direction of spin current, 
the unit vectors X and Y are perpendicular to Z, and J denotes 
the spin current amplitude vector, which is an axial vector par- 
allel to the spin polarization direction. Now the z and Z axes 
form a special plane. If the system has reflection symmetry 
with respect to this plane (e.g., the system is spherically sym- 
metric or the plane is along a special crystal direction of the 
semiconductor), the symmetry properties of the relevant quan- 
tities under reflection with respect to the z-Z plane will impose 
further constraint on the interaction and significantly simplify 
the Hamiltonian. Under the reflection, the relevant quantities 
transform as 





q 


h 


Jll 




Reflection with z-Z plane 











(12) 



where Jy is the component of J in the plane and is the 
perpendicular component. By the table above, it is evident 
that to keep the efl'ective Hamiltonian invariant only the in- 
plane component of Jy would couple with the ql^. Without 
loss of generality, let Y be perpendicular to the z-Z plane, the 
efl'ective Hamiltonian reads 

=A,qI,Jz+A2qhJx, (13) 

in which two coupling constants Ai and A2 are to be deter- 
mined by microscopic calculation. Alternatively, the Hamil- 
tonian can be expressed in a form independent of the choice 
of the X and Y axes as 

9Y^^^ = fi^/,/z + 6^/zZ-J-z, (14) 

which is the same as Eq. ([9]), but does not require the spherical 
symmetry of materials. 

The physical efl'ect of the efl'ective coupling can be ex- 
tracted from the linear optical susceptibility, 

X^,v +Xl,, = (l/eo)d^9{,ff/(dF;dF,), (15) 

where 60 is the vacuum permittivity. Thus we get an opposite 
susceptibility for opposite circular polarization in presence of 
a spin polarization or a pure spin current 

X?^ = -x^l = (l/46o)foz • S, (16a) 
/!l = -X^ll = (q/46o) (fi/z + f2Z • J • z) . (16b) 
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The effective energy shift resulting from the Hght-spin or 
Hght-spin current interaction means a phase shift in the Hght 



observed in the far-field. Eq. (16a) is nothing but the con- 
ventional Faraday rotation in magnetooptics,^^ Eq. (16b) in- 
dicates that a pure spin current would produce a circular bire- 
fringence eff'ect. This new eff'ect of a pure spin current may be 
dubbed "spin current Faraday effecf ^^ because of its similar- 
ity to the conventional Faraday rotation due to magnetization, 
with awareness that a pure spin current, however, bears no net 
magnetization. 



B. Second-order nonlinear optical effects 



or X must be zero, so 



xf^ =Jx (xiXXY -h JC2XYX -h JC3YXX -h ziZZY 

-hZ2ZYZ -h Z3 YZZ -h jYYY) , (20) 

with seven independent parameters to be determined. Similar 
symmetry analysis can be applied to JyYZ. Such unique po- 
larization dependence of the second-order optical susceptibil- 
ity can be used to distinguish the longitudinal and transverse 
components of a spin current, and also to single out the spin- 
current signature f rom t he eff'ects of the material background 
or a charge current.l^^l^ 



The second-order nonlinear optical eff'ect such as sum- 
frequency process is characterized by a second-order nonlin- 
ear susceptibility x^^"* via 



(17) 



where Fi and F2 are the two optical fields, P is the induced 
polarization, and x^^"* is a rank- 3 tensor. Under P operation, 
Fi , F2 and P reverse the sign, which means x^^^ is zero if 
the system has ^-symmetry. A pure spin current breaks the 
!P- symmetry, results in a nonzero and makes the second- 
order nonlinear optical process possible. 

In general, as a rank-3 tensor has 27 independent com- 
ponents 

= XxxxXXX ^XxxyXXY + . . . ^XzzzZZZ. (18) 

But the symmetry properties of the spin current and the sys- 
tem will impose constraints on x^'^\ reducing the number 
of independent parameters For a longitudinal spin current 
/zZZ, in which the spin polarization is parallel or antipar- 
allel to the current direction, the spin current is reversed 
under the reflection with respect to the X-Z plane, so that 
XxxxXXX +xxxyXX(-Y) + . . . + XzzzZZZ = -x^^^ and the 
terms with direction Y appeared even times (twice or zero 
times) must vanish. Similarly, the longitudinal spin current 
is reversed under the reflection with respect to the Y-Z and X- 
Y planes, so the terms with direction X or Z that appear even 
times must be zero. Moreover, the longitudinal spin current 
should be invariant under the tt/ 2 rotation with respect to its 
current direction, therefore xxyz = -Xyxz, Xyzx = -Xxzy, 
XzxY = -XzYX- With all these constraints, the sum- frequency 
susceptibility induced by a longitudinal spin current JzZZ can 
be expressed as 

xf^ =Jz [0^1 (XYZ - YXZ) -h Qf2(YZX - XZY) 

-^asiZXY - ZYX)] , (19) 

where there are only three independent parameters at (i = 
1, 2, 3). For a transverse spin current JxXZ, in which the spin 
polarization is perpendicular to the current direction, the cur- 
rent is reversed under reflection with respect to the X-Z plane, 
but invariant under refection with respect to X-Y or Y-Z plane. 
Then the terms that contain even times of Y or odd times of Z 



III. MICROSCOPIC CALCULATION 

To quantitatively determine the linear and the second- 
order nonlinear optical efl'ects of a spin current, we will per- 
form the microscopic calculation for a pure spin current in 
a bulk direct-gap semiconductor using the standard perturba- 
tion theory.^^ We employ the eight-band model.^^ We as- 
sume that the pure spin current result from a non-equilibrium 
distribution of electrons in the conduction band. Namely, as 
shown in Fig. [TJa), a small portion of non-equilibrium elec- 
trons with opposite velocities near the Fermi surface have op- 
posite spin polarizations, which is similar to the situation in 
Ref . 8 The optical interaction includes the interband transi- 
tions and the intraband acceleration of electrons and holes. To 
avoid real absorption of light, the light frequencies are chosen 
to be below the absorption edge in linear optical efl'ect, and 
the sum frequency is below the two-photon absorption edge 
in the second-order nonlinear optical eff'ect. 



A. Model 

We consider an n-doped direct-gap semiconductor of GaAs 
as a model material. Since other bands are separated far 
away in energy, we assume the near-gap optical interactions 
in GaAs involve mostly the eight bands around the fundamen- 
tal gap, including the conduction band (CB), the heavy-hole 
(HH) band, the light-hole (LH) band, and the spin-orbit split- 
ofl' (SO) band, each of 2-fold degeneracy [Fig.[TJa)]. 

Near the F-point of the Brillouin zone, the energy dis- 
persion of the CB and the SO electron is almost parabolic 
and isotropic, which can be respectively written as Eep = 
p^/(2me) + Eg and Etp = p^/(2mt) + Eso, where nte/t is 
the efl'ective mass of the CB/SO band. Eg is the fundamental 
band gap with the energy zero at the top of the valence band, 
and Eso is the split-off' energy due to the spin-orbit coupling. 
Hereafter the Planck constant fi is taken as unity. For the HH 
and LH bands, since the energy band anisotropy will not aff'ect 
the symmetry analysis as shown in Sec.|ll|( The 4-fold rotation 
symmetry will be retained even considering the anisotropic 
dispersion), we will neglect the anisotropy eff'ect in this Sec- 



tion, and take it into account separately in Sec. IV B Thus we 



express the isotropic Luttinger-Kohn Hamiltonian H[j^ for the 
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FIG. 1 : (color online) (a) Schematic band structure of the eight-band 
model near the T point of an n-doped III-V compound semiconduc- 
tor, and illustration of the spin non-equilibrium distribution of elec- 
trons for a pure spin current, (b) & (c) Selection rules and relative 
dipole moments from the HH and LH, SO bands to the CB. 



HH and LH bands near the band edge ai 



m 



1 

2mo 



7i + 



- 272 (V • K)^ 



(21) 



where K is a spin- 3/2 for the total angular momentum of an 
electron in the HH and LH bands, yi and 72 are the Luttinger 
parameters, and mo is the free electron mass. The isotropic 
Luttinger-Kohn Hamiltonian can be diagonalized with the 
spin-3/2 quantized along the direction of the momentum p. 
The HH band with magnetic quantum numbers ±3/2 has the 
energy dispersion Ehp = (71 - 272)/?^ /2mo = p^/2mh', and the 
LH band has magnetic quantum numbers +1/2 with the dis- 
persion relation as Ei^ = (71 + 2y2)p^/2mo = p^/2mi, where 
ruh/i is the effective mass of the HH/LH band. If the HH-LH 
splitting is neglected further, the HH and LH bands become 
a 4-fold degenerate spin-3/2 band and the spin quantization 
direction can be chosen independent of the momentum. 

The Bloch state of a CB electron with momentum p is 
|j/^^(p)) = ^1 p|0) = ^^P'^|+)p, where |0) represents the vacuum 

state and ^ denotes an creation operator that produces an 
electron in CB with spin +1/2 and momentum p. The Bloch 
state of an electron in the valence band with momentum p 
is |C(P)) = Vl^;p|0) = e'nj^m}^, in which j = 3/2 and 
m = +3/2 stands for the HH band (a = h), j = 3/2 and 
m = ±1/2 for the LH band (a = /), j = 1/2 and m = +1/2 for 
the SO band (a = t), and V;,m;p denotes the annihilation oper- 
ator for an electron in the corresponding valence band. Then 
the non-interacting Hamiltonian is 

+^rp?IpW). (22) 

^ f/t 



where the hole operators are defined as /z+_p = V^i2 + 



t 

3/2,±l/2;p 



and f+ _p = V, 



t 

■ l/2,±l/2;p 



3/2,+3/2;p' 

. It should be pointed 



out that here the angular momentum K is quantized along p so 
that the spin-orbit coupling in the valence bands has already 
been included. 

The initial state of the system (before optical excitation) 
is characterized by a density matrix po. We assume that the 
system has translation symmetry and initially there is no hole 
in the system, so we have 



Tr[po^^,^^yk'] = ^Kk'ffivM- 



(23a) 
(23b) 



The spin current, which results from the non-equilibrium dis- 
tribution of CB electrons, is expressed by Eq. ([3]), where the 
velocity and the spin polarization of an electron with mo- 
mentum p is respectively given by Vp 
(1/2) Y,j2v ^fivffiv,p with a denoting the Pauli matrices. 



Vp^gp and Sp 



B. Linear optical effects 

The direct interaction between a light beam and a semi- 
conductor is the dipole interband optical transitions. Only 
through the spin-orbit coupling in valence bands, may the 
light beam interact with the spin of electrons. 

For the dipole interband transition [Fig.[TJb)&(c)], the po- 
larization density operator reads^^ 

PinterCl") — ~~ ^cv j l^//,p^//,-p^/i,k "I ^llyu,p^,-p^/i,k 

J J - - /2 . . 
+ ^zp?^,_p^^,kj ^^P"-^^" + H.c, (24) 

where n+,p = + (xp + iy^^ / V2 denotes the right/left circular 
polarization about the momentum direction p, Zp = p//?, and 
fl = -[I. As will be discussed in Sec, |IVB[ the momentum 
dependence of the dipole moment has no significant eff'ect, 
so here we assume the interband dipole moment Jcv indepen- 
dent of the momentum. With the dipole interaction with a 
fight Hiit) = - J Pinter(r) ' F(r, 0<ir, the light-matter interac- 
tion Hamiltonian in the rotating wave approximation can be 
explicitly expressed as 



Hi = exp ^ iojqte^^^^e^^^ 
V//,k 



^i(Oexp 



^cv Yj ^y^y ' (>^//,p^//,-p^//,q+p + r-'^yU,p(u,-p^/i,q+p 



yU,V,p 

f2 



Zp///,-p^//,q+p 



Il//,p ^//,-p ^yu,q+p 



"I" ^^^p^yU,-p^A/,q+p I H.C 



(25) 



Under the condition that the optical interaction strength is 
much smaller than the detuning of the light from the valence 
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band to the Fermi level (the perturbation regime), the effec- 
tive energy due to the dipole interaction can be derived by the 
second-order perturbation as 



Hff = Tr 



(26) 



Such effective coupling between a spin current and a polar- 
ized light beam on the one hand can be regarded as the fre- 
quency shift of the light in the presence of the spin current, 
and on the other hand can be considered as the energy change 
of the semiconductor system under the driving of light beam. 
The second-order perturbation means that there are two virtual 
optical transitions induced by the electric field of the light: 
one creating an electron-hole pair and one annihilating the 
electron-hole pair. The virtual excitations cause no real op- 
tical absorption but a phase shift, indicating that the effective 
coupling is real. The optical effect of the spin-current can be 
understood as the Pauli blocking in the transition involving 
different spin states. With this picture in mind, the following 
microscopic calculation, though lengthy, is quite transparent. 




FIG. 2: (color online) Physical picture for the microscopic mecha- 
nism of the spin current Faraday effect. The virtual transition energy 
is the same for +Sk when neglecting light momentum q. 



1. Physical Picture 

The physical picture for the microscopic mechanism of the 
spin current Faraday effect is rooted in the fact that a spin will 
induce a Faraday rotation like a magnet. In Faraday rotation, 
a linearly polarized optical field F induces a polarization as a 
rotation about the spin. 



p(i) 



Fxsk 



(27) 



where Sk is the spin polarization associated with the state of 
k, £'k is the resonant optical transition energy. This natu- 
rally explains Faraday rotation due to spin polarization as in 



Eq. ([36a). 

For a pure spin current, we first consider only a pair of 
spins, Sk at momentum k, and -Sk at momentum -k in the 
CB [Fig. [21. This pair can be viewed as a generator of pure 
spin current. Sk gives rise to a Faraday rotation of Pj^ ^ oc 
F X Sk/(6L) - £^+k); while -Sk leads to a Faraday rotation of 
P^k ^ F X Sk/(6i; - £^-k). Therefore, the Faraday rotations 
caused by the pair of spins cancel each other in the vertical 
optical transition. However, when the effect of the small light- 
momentum is taken into consideration, the excitation energy 
at ±k will shift respectively to E+\^ £'q+k ~ £'±k ± q • Vk 
[Fig.[2|, and ±Sk will induce different Faraday rotations due 
to opposite energy variation. Up to the first order of q, the po- 
larization is P^^^ oc FxskVk-q/(^^-£'k)^, where ^SkVk is just the 
spin current tensor contributed by the pair of electrons. This 
explains the ^-dependence of the spin current Faraday effect. 
More generally, the hole state wavefunction is also changed 
when considering the light momentum, which causes extra 



2. Effective Hamiltonian by SO-CB transitions 

To better understand the microscopic mechanism of the 
light-spin current coupling, let us first derive the effective 
Hamiltonian contributed solely by the transitions between the 
CB and SO bands. The SO band electrons has 2-fold degen- 
eracy, and the spin states as well as the selection rules for the 
interband transitions, like the CB electrons, are independent 
of the momentum [Fig.[TJc)]. 

We first consider a single electron with momentum k and 
spin polarization Sk. The spin current contributed by this elec- 
tron is Jk = ^SkVk with the velocity Vk = k/nie. It is conve- 
nient to define the spin basis states along the spin polarization 
direction. In such chosen basis, the spin density matrix of the 
electron is diagonal. With the population in the spin-up and 
spin-down states denoted as /+ and /_, respectively, the spin 
polarization is ^k = (/+ - /-)/2. The interband transitions 
|l/2, ±l/2)k |=F)k couple to a field with circular polariza- 
tion =F(ei ± /e2)/ V3, and the vertical inter-band transitions 
|l/2, ±l/2)k |±)k couple to a field of linear polarization 
63 / V3 , where the coordinate system is so defined that 63 is 
along the spin polarization direction of the electron consid- 
ered. Summing up all possible inter-band transitions, the en- 
ergy shift of this electron due to coupling to an optical field F 
is 



^80 _ _ 1 , . ,2 V (l-/±)F*-(ei + /e2)(ei + /e2)*-F 



(l-/JF*-e3e;-F 

- Et-k+q - Ee^ 



(28) 



Berry phase effects [terms proportional to 1 lEp in Eq. ( 37 )]. 



where the factor (1 - /±) accounts for the Pauli blocking of 
the interband transitions. The second term, which is related to 
vertical transitions caused by a linearly polarized field, does 
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not depend on the spin polarization, so it can be dropped as 
the background. With expansion to the first order of q and 
omission of the background terms, the energy shift becomes 



2 

+ -i\dc^ 



|2 >ykF* -(6162- 6261) -F 
2 ^kF* -(6162- 6261) -Fq-k 



(29) 



Since (6162 - 6261) • F = F x (ei x 62) = F x 63, the physical 
meaning of this coupling is transparent: the linear-polarized 
optical field will tilt about the spin, which is essentially the 
Faraday rotation with spin playing the role of a magnet. The 
summation over the momentum space gives 



3 At 3 erritAf 



qI,Z']'Z, (30) 



where is the light detuning from SO band to the Fermi sur- 
face. 



3. Effective coupling by transitions between HH/LH and CB 

If the HH bands and LH bands are assumed degenerate, 
the quantization direction of the 3 /2-spin of the HH and LH 
can be chosen arbitrarily and the eff'ective Hamiltonian are ob- 
tained in a similar way to that contributed by the SO-CB tran- 
sition. However, with the HH-LH splitting considered, the 
quantization direction of the hole states depends on its mo- 
mentum, thus, with the trivial background omitted, the eff'ec- 



tive Hamiltonian can be derived explicitly as 
"Hf^ =l^cvP [h (XX - yy) -h ly (xy -h yx) -h /, (/xy - /yx)] : 



if 

D Hi, 

V2 



n/i,pn^,p + 2zpz; 



3 ffipfip,q+p ^ 



eq+p 



+ E^ 



/-p 



(31a) 

(31b) 
(31c) 



where yUp indicates the spin moment quantized along p. The 
terms in the equation above stand for diff'erent physical pro- 
cesses as follow. Term ( [31a| ) accounts for the HH-CB tran- 
sitions where an (virtually) absorbed photon will be emitted 
with the same circular polarization conserving the electron 
spin. The LH-CB optical transition can be either circularly po- 
larized or linearly polarized. In ( |31b| ), the absorbed and emit- 
ted photons in the virtual LH-CB transitions have the same 
polarization with the electron spin conserved. In ( |31c| ), the 
optical polarizations involved in the LH-CB absorption and 
emission are changed, leading to an angular momentum trans- 
fer between the light and the electron, while the total angular 
momentum is still conserved. 

The spin-independent population terms such as /++ p p 
are related to a change in the background refraction index, but 
not to the spin polarization and spin current, we will drop the 
spin-independent population terms and keep the spin polariza- 
tion terms Sp only in the eff'ective coupling as 



^eff"^ = l^cvl' [h (XX - yy) + /, (xy + yx) + //, (xy - yx) ] : / ^ 

p 

^P'P+q(yp^p ~ Zpyp) ^p,p+q(Zp^p ~ XpZp) 



-e,p+q 



+ E, 



■ flOJa 



/zp,p+q(xpyp ypXp) 1 



-e,q+P 



+ E, 



(32) 



where /e.,p = Sp-e/. Using the anti- symmetric tensor = t/yy^e/e^ey^ which is invariant under orthogonal coordinate transformation. 



we can express Xpyp - ypXp = Zp , ypZp - Zpyp 



Xp • (o and ZpXp 



yp • , whereby the terms associated with the 



electron spin polarization form an anti- symmetric tensors. Noticing that the contraction between the anti- symmetric and the 
symmetric tensors associated with and ly must vanish, and also that the eff'ective Hamiltonian must be real, we have 



^eff'^ = -Kvl'/z(xy-yx):^ 

p 



Sp+q • ZpZp • S' 



Sp+q • ZpZp • (o 



2 

-h -- 



-h Eir, - flCOa 3 £gq+p -h £/p - flCD, 



=2\d,,\'I, 



E 

P 



■eq+p 



(Eh 



3 Eeq+p + Eip — flCDq 



(33a) 
(33b) 



By expanding to the first order of q, we have £'eq+p ~ Ee^ + q • Vp£'gp, and Sp+q Sp q • VpSp. By using VpZp = Vp(p//7) = 
P^lp - pp//?^ = (XpXp -h ypyp)/;?, and Vp(ZpZp) = (XpZpXp -h ypZpyp -h XpXpZp -h ypypZp)//?, we obtain the eff'ective Hamiltonian 
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as 



2 ^ 



/ Me 



Me 



,2 Me 1 



3 m/A^ 2A/£^ 2A;,^^ 



(34) 



r 



where Ep is the Fermi energy, A/^// is the Hght detuning from 
the HH/LH band to the Fermi level, respectively [see Fig. 
1(a)]. The first term in Eq. ([34]) results from the spin po- 



larization, while the other terms result from a spin current. 
When neglecting the HH-LH splitting by letting A/^ = A/ and 



Mh = mi, Eq. p4\ is reduced to a expression similar to Eq. ( |3Q| ) 
but with a minus sign(m^ ruh, A^ Ah). This reduction 
confirms that the eff'ective Hamiltonian from the HH/LH-CB 
transitions can be derived as easy as that from the SO-CB tran- 
sitions if the spin quantization direction in HH/LH band can 
be chosen arbitrarily. Moreover, if there were no spin-orbit 
coupling in the valence bands, i.e., the HH, LH, and SO bands 
had the same eff'ective mass and the same band-edge energy, 
the coupling between a spin current and a light would vanish. 

Finally, once the spin distribution is specifically given, the 
total eff'ective Hamiltonian will be determined. We assume 
that the electron spin distribution around Fermi wavevector 
deviate only slightly from the equilibrium distribution. More 
specifically, we suppose the spin distribution has the form 



Sp = No + Ni/(/7)cos^p, 



(35) 



where is the angle between the momentum p and the cur- 
rent direction Z. Such a distribution is the usual case for weak 
currents. A straightforward integration over the momentum 
space gives [see Appendix [Aj 

^l^ = ^oh^-S, (36a) 



with the coupling constants 
_ 2 

"a; 

2me 



^2 = 



+ 



e \5ot/,A^ 5ot/A^ 5ds.hEF 5AiEf 
p/ Anie Smg 4me 1 



\dc 



(37a) 
(37b) 

1 



e \5mhAl ISm/A^ 3mtAj SAhEp SAiEp )' 

(37c) 



For a spin distribution diff'erent from Eq. ( [35] ), as can be seen 
in Sec. IV D, the coupling constants shown above will only be 
changed quantitatively, which further confirms the symmetry 
analysis in Sec. II [see Eq. ([9])]. 



C. Second-order nonlinear optical effects 

The linear optical eff'ect of spin currents is weak since the 
photon current involves the small light momentum. If we re- 



place the light momentum by another optical field, the cou- 
pling can be greatly enhanced by a factor of F2 • v/q • v. As 



shown in Sec. II B the second-order nonlinear optical eff'ects 
of spin currents is rooted in their unique physical nature and 
spatial inversion- symmetry breaking. Specially, noticing that 
a longitudinal spin current, is a chiral quantity, we envisaged 
that it could be probed by the chiral sum-frequency optical 
(SFG) spectroscopy which was recently developed to detect 
molecular chirality.^^ If otherwise measured in linear op- 
tics, the eff'ect of the chirality relies on the small magnetic 
moment of the molecules, and in turn on the small wave vector 
of the probe light, similar to the case of linear optical eff'ects 
of spin currents. 



1. Physical picture 

The nonlinear coupling between a spin current and light 
contains three processes: one virtual interband transition cre- 
ating an electron-hole pair, one intraband transition acceler- 
ating the electron or the hole, and one virtual transition in- 
ducing the combination. The physical picture for the micro- 
scopic mechanism of the second-order nonlinear optical ef- 
fects of spin currents is similar to the linear optical eff'ect. A 
spin will induce a Faraday rotation P^^^ oc F x - E\^). 




FIG. 3: (color online) Physical picture for the microscopic mech- 
anism of the second-order nonlinear optical effects of a pure spin 
current. The second light F2 will accelerate the electrons (or holes). 
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The Faraday rotations due to the pair of spins of Sk (at mo- 
mentum k) and -Sk (at momentum -k) cancel each other in 
the vertical optical transition. Instead of considering the small 
light-momentum in the linear optical effect, we add another 
optical field F2. The spin will experience an intraband ac- 
celeration by this optical field and the transition energy will 
be changed to ^±k ^ ^±k ± / e\i, • ¥26-'"^'^ dt [Fig.jsj. The 
physical meaning of ^Vk • F2 is clear that is the power aone by 
the field to the electron. Therefore, ±Sk will induce diff'erent 
Faraday rotation due to opposite energy modification 

P(2> cx Fi X SkCVk • ¥2/[(aJi +CD2- E^XoJi - E^)aj2l (38) 



This gives the second-order nonlinear optical eff'ects of spin 
currents. 



2. Microscopic calculation 

The second-order nonlinear susceptibility can be obtained 
straightforwardly through the standard perturbation method as 
shown below. Here we take the SFG as an example of the 
second-order nonlinear optical eff'ects of spin currents. 

The dipole density operator for the intraband transition 
reads'^ ^ 



Pintra(r) ^ ^ ^^.,p^//,k<A^%/^)k + ^ V),^..p Vj,m;k</, ^'Ip J, w)k 
k,p ^fi,fi'^± f,m';j,m 



y^^/pr-/kr 



(39) 



With the input optical field consisting of several frequency components F(r, t) = Yjj^i,2 ^j^ '^^^ + c.c., the light-matter interaction 
Hamiltonian is Hi(t) = - J P(r) • F(r, t)dr, where P(r) = P(r)inter + P(r)intra- Explicitly, we can write 



= - (d -h d) • I ^ Fje-'"'^' -h c.c. j, 
^ 7-1,2 / 



(40) 



with 



6 = - <v ^ (n^,k/i/i ,-k^/i,k + (1/ V3)n^,k//i ,-k^/i,k - V273zk//i -k^/i,k - V273n^,k?/i ,-k^/i,k + (m/ V3)zk?;i -k^/i,k) , (41a) 



d=/^^Q]^i,p^M + Z^I-;p^^'^^^ Z ^i',k^M<y^'lkVkl//)k + ^ V]^,.kV7,m;kO''^'lkVkl7,m)k I (41b) 

k,p ^ {!-+ j,m ' k ^{i,ii'-± j,m,m' ' 



denoting the inter- and intra-band polarization operators, re- 
spectively. D and are the positive- and negative-frequency 
components of the inter-band polarization operator, respec- 
tively. The first part of the intra-band polarization is the usual 
acceleration term. The second part, which has the form of 
non-Abelian Berry connections (similar to vector potentials), 
accounts for the variation of the spin quantization direction 
with acceleration of an electron. It is necessary to include the 
Berry connection term for the gauge-invariance of the intra- 
band polarization. The explicit form of the Berry connection 
term depends on the choice of the local coordinate (Xp, yp, Zp) 
at momentum p. In Appendix [B] we present an example for 
the Berry connection in a specific convention. 

We adopt the interaction picture for calculating the SFG. 
The second-order polarization response obtained by the stan- 



dard perturbation theory is 

F^\t) = - f df f j/Tr[D(O[A(^0,[A(^"),Po]]], 

where D(0, Hi are operators in the interaction picture. We 
consider the case that (1) the sum frequency oj = oji -\- 0J2 is 
near resonant with the band-edge, so the positive-frequency 
component D(0 dominates the optical process; (2) the intra- 
band dipole moment must be considered for the contribution 
by the spin current; and (3) no holes exist in the initial system, 
so the inter-band excitation has to be involved (caused by D^). 
With all these considerations taken into account, the second- 
order response of interest is 



p(2)(^) = - J J df' e-'""''' Tr{f)(t)¥2 • f)\f) Fi • d(^'0,Po|) 



•D^(r")Po]) 



+ {Fi,a;i F2,a>2 
+ {Fi,a;i <-> F2,a;2}. 



(42a) 
(42b) 
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The physical meaning of Eq. ( 42 ) is clear: Eq. ( 42a ) corre- 
sponds to the driving of the electron population (at fol- 
lowed by inter-band excitation (at f) and emission (at t); 
Eq. ( |42b| ) corresponds to the process in which an electron- 
hole pair (created at is driven by an external field (at f) till 
its emission (at 0- 

When the HH-LH splitting is neglected, we have a sim- 
ple microscopic calculation as discussed in Ref. 21 , in which 
the spin quantization for valence band states and the selection 
rule for interband transitions are independent of its momen- 
tum. Beyond such an approximation, the calculation of P^^^ 



through Eq. (42) is lengthy, but only quantitatively modifies 
the results. So we will only list the result in the Appendix [Ej 
and the details are shown in the Supplementary Information. 




IV. DISCUSSIONS AND NUMERICS 



A. Faraday rotation of a spin current and spin polarization 



FIG. 4: (color online) The geometry for measuring a spin current, in 
which the spin current is along Z-direction and the red arrow denotes 
the propagation direction of the light beam. 



The Faraday rotation angle is expressed as 

= OJql(x++ -X—)l (4^<^) . 



(43) 



where / is the light propagation distance, n is the material re- 
fractive index, and c is the light velocity in vacuum [Appendix 

Pure spin current. For a spin current configuration as 
shown in Fig. |4] where a light comes in with a zenith angle 
P and an azimuth angle y, the Faraday rotation angle due to 
diff'erent components of JZ is 



^^.Vx) = 4''^^f2SinySC0Sr, 

ef^Jy) = -S^^^Y^i sin^Psiny COS y(n^ 



(44a) 

sin^yS)-^/^ (44b) 

sinV)"^^^ 
(44c) 



where = n^l/2n€oA^. The dependence of the rotation an- 
gle on the incident angles for Jz, Jy and Jx components of a 
pure spin current are shown in turn in Fig. [5] (a), (b) and (c). 

Net spin polarization. The net spin polarization also causes 
the Faraday rotation. With the incident light of zenith angle 
the Faraday rotation angle equals 



C(S) : 



(InllU^nX) (^oz • S) . 



(45) 



Spin polarization has both the normal and parallel components 
with respect to the sample surface S = + Sy . For the normal 
component S^, the rotation is independent of yS, 



^f(S^) = 7r^o^^L/46o^^, 



while for parallel component Sy 



ef(S\\) = (7i^oSiiL/4eoA) smficosy(n^ - sin^ fi) 



-1/2 



(46) 



(47) 



In general, the angle dependence of Faraday rotation can 
be used to distinguish a pure spin current from a spin polar- 
ization. However, in many materials ^ » 1 > sinyS, both 




FIG. 5: (a)-(c) The Faraday rotation amplitude of spin current com- 
ponents Jz, Jy and Jx as functions of the incident angles of the 
light beam, (d) The dependence of (5p ^ on the light wavelength A. 
Parameters are chosen similar to those in Ref. 8 Eg =1519 meV, 
E^Q =341 meV, the doping concentration is 3 x 10^^ cm~^, the ef- 
fective mass (in units of free electron mass) of the HH, LH, SO, 
and conduction bands is in turn 0.45, 0.082, 0.15, and 0.067, the 
dipole dcv - 6.7 eA, n - 3.0, L - 2.0 jim, Ep - 5.3 meV, and 
Jx^Jy^Jz^ 20 nAyL/m-2. 



\S||) and 0^^\jx) have nearly the same angle dependence, 
which is proportional to sin /S cos y. As there is inversion sym- 
metry diff'erence between a pure spin current (P = odd) and 
a spin polarization (!P = -h, even), a pure spin current would 
have a sign flip at reflection while a spin polarization would 
not. Therefore, the Faraday rotation angle of a pure spin cur- 
rent vanishes through reflection, while the rotation angle of a 
spin polarization will be doubled. This difl'erence can be used 
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to distinguish the effect of a spin current from that of spin po- 
larization. 

For the reaUstic case in Ref . |8l the vanishing Faraday sig- 
nal is reported in the middle region where the spin current 
flows without net spin polarization. We explain it with the 
fact that in the experiment Z • z = and /z = 0.1^ With the 
experimental configuration shown in Fig.|4j the rotation angle 
6^^^(j3,y) oc sin cosy. The maximum Faraday rotation 
angle is reached when p ^ njl and y ^ 0. The dependence 
of maximum Faraday rotation angle on the light wavelength 
is plotted in Fig. [5|d). For the specific example shown in 
Fig.|5jd) with light wavelength around 800 nm , the maximum 
Faraday rotation angle is 0.38 yurad. Such a Faraday rotation 
angle, though still small, is measurable in experiments. 



The calculation is lengthy. Here we omit the details but just 
give the terms with i = 3, = /u' = -\- and / = -3, jj = /j' = - 
explicitly, which is proportional to 

(\Al^\\x^ + lA^ jVpyp + j2zpZp)(A^,p + /__,p) (53a) 
-h 2/[3(A^^A2p(xpyp - ypXp) -h 2/3(A2^A^p(ypZp - Zpyp) 
+ 2/3(A^^A^;)(ZpXp - XpZp)]^. (53b) 



The term ( |53a| ) is just a background. The term ( |53b| ) is the 
total anisotropy-tensor, which couples to 1^ only. This result 
confirms the symmetry analysis in Sec. [Ill 



B. Effects of valence band anisotropy 



C. Second-order nonlinear optical effects 



In the derivation above, we have neglected the anisotropy of 
the valence bands. Now we examine the eff'ect of the valence 
bands anisotropy. The anisotropic valence band Hamiltonian 
takes the form 



The independent parameters of the susceptibility of spin 
current in a bulk GaAs in Eqs. ([19]) and ([20| are listed in 



HtK =7^ \(n + 572/2) - 273 (V . Kf 

+2(r3-r2)(vX' + c.p.)], 



(48) 



where the (73 - 72) term describes the anisotropy. The 
anisotropy is usually small. The eigenfuctions of //^^ are 

\i/ft}= Yj ^/|3/2,j), / = ±3,±1, (49) 

7-±3/2,±l/2 

where the basis states |3/2, ±3/2) and |3/2, ±1/2) are explic- 
itly given in Appendix pi and are coefficients satisfying 
V^a = Of*, with U = -icTx ® o-y. The eigenstates 10^+3) and 
|j/r+i) have eigenvalues Eh and £/, respectively. The dipole 
density operator can be explicitly written as 

P(r) = -eY,e-'''-'^'''-%,-^e,,M 

+ /+,_p^^,k<<A-i|r|A/) + /i+,-p^;,,k(<A-3|r|yu)] + H.c, (50) 

where = |±) denotes the CB electron state with spin ±1/2, 
and the operators h+ and /+ annihilate {1/^+3} and |</^±i), respec- 
tively. By using the fact p/mo = dr/dt = (vHq - Hor)/i, we 
get 



(51) 



J^x,y,z 



where M+3/+1 p = -ie/mo(Ee^p - Ehu^). The detailed expres- 
sion for A^^ can be found in Appendix C. The eff'ective Hamil- 
tonian then reads 

I j,j'=x,y,z 
^ cr,cr' ,p,i,fi,fi' 

iMrJ A. A. Jul n- 

I ''P| l,IJ. ^yU'JP*' P J 



(52) 



Appendix |E| For the sake of simplicity, we neglected the 
anisotropy of the valence bands. We also neglected the 
Coulomb interaction, since it is largely screened in the n- 
doped material. These approximations, according to the sym- 
metry analysis, would only quantitatively modify the results. 
The bulk inversion asymmetry would cause a background 
second-order susceptibility, which is indeed strong but can 
be well separated from the spin-current eff'ect by ac modu- 
lation of the current and the phase-locking detection tech- 
nique. Two representative results of the calculated suscep- 
tibility spectra are shown in Fig. |6] The other terms of the 
susceptibility tensor (not shown) have similar frequency de- 
pendence and comparable amplitudes. As a specific exam- 
ple, a transverse spin current 20nA/yum^ has a susceptibility 
-^f^^ ^ 4.8 X 10-i2^su (or 0.2 X lO'^^ -^^ 3J ^j^j^^) 
for input frequencies 6i;i = 100 meV and 0)2 = 1400 meV or 
0.25 X 10"^^ esu for oji = 0J2 = 750 meV (corresponding to 
the second harmonics generation). 

The SFG of spin current can be straightforwardly extended 
to other second-order optical spectroscopy such as diff'erence- 
frequency and three- wave mixing. 



V. CONCLUSIONS 

In summary, with the systematic symmetry analysis in gen- 
eral and the microscopic calculation under realistic condi- 
tions, we have shown that a pure spin current has a measur- 
able circular birefringence eff'ect and a sizable sum- frequency 
susceptibility. With universality of the method guaranteed by 
the symmetry principle and without requirements of special 
structure design and fabrication, the linear and nonlinear op- 
tical spectroscopy can be applied to study a wide range of 
spin-related quantum phenomena such as in topological insu- 
latorJ^^Jtm A wealth of physics connecting spins and photons 
and technologies synthesizing spintronics and photonics may 
be explored. 



12 

























J 




i 












j 



1200 



750 1000 o>v 




1 200 0^^^ 



. onn 1100 

750 1000 



And the angle average 



^ ' ZpZpZpZp 



4;r/ 



^ * ZpZpZpZp^O 



+cos2 sin^ ^p cos2 0p (XXZ + XZX + ZXX) 



+ cos2 ^p sin^ ^p sin^ 0p (YYZ + YZY + ZYY) 



+ cos4 ^pZZZ 



^ [P^Z + XZX + YZY + ZZZ + ZP^) . (A3) 



For a spin distribution of Eq. ( [35] ), the total spin polarization 
and the spin current is respectively as 



P P 

J = y jp = — y Ni/(/7)/7z . zpzp 

p p 

= — ^Ni/(/?)/7Z-z^ 
^ p 
NiZ e 



3 



(A4a) 



(A4b) 



FIG. 6: Representative results of the sum frequency susceptibil- 
ity, (a) -Xyzx ^ longitudinal spin current, and (b) -Xyyy 
due to a transverse spin current, as functions of the optical frequen- 
cies. Parameters are chosen similar to those in Ref. 8 (same as in 
Fig. [5]). The dielectric constant 6^ = 10.6, and the spin current 
Jx = Jz = 20 nA/yum^. 
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Appendix A: Coordinate basis 

We choose a global coordinate system (X, Y, Z) and define 
the local coordinates as 

§p =Xp = X cos ^p cos 0p + Y cos ^p sin (pp-Z sin 6^, (Ala) 
0p =yp = -X sin 0p + Y cos 0p, (Alb) 
p =Zp = X sin ^p cos 0p + Y sin ^p sin 0p + Z cos ^p. (Ale) 



Also, we have 



The angle average of the tensor 



(A2) 



^ ZpZp • {Nif(p)pZ • ZpZp) 
p 

= — ^ /(/?)/?Z • ZpZpZpZp • Ni 



e 

Me 

e 

Me 



1 e 

3 Me 



/^2)z + XZX + YZY + ZZZ + Z/(2) 
X • Ni 



(A5) 



Appendix B: Berry connection 

The band edge state of CB are 

l+/->p = l^>^IT/i>p, 



(Bl) 



with 15" ) being a periodic s-wave orbital wavefunction which 
is isotropic in a unit cell, and | T / >l)p denoting the spin eigen 
state parallel/anti-parallel to the momentum. 
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Similarly, the band edge states of the valence bands are 
l^>p + i\Y\ 



3 3\ 
3 1^ 

r 21 



V2 



' I T>p, 



2 |X)p + i\Y)„ 



|Z>P ® U>p + 



\X}p - i\Y\ 



2'"^2/ 

1 _1^ 

2' 21 



V3 



V2 

|Z)p®IT>p 



V6 

\x)p - /|r)p 

V6 

i>p, 

|X)p + /IDp 
V3 



®IT) 



1 , \X)p-i\Y)p 

.^|z>p«U)p ^ 



(B2a) 
(B2b) 
(B2c) 
(B2d) 
^ I i)p, (B2e) 
^IT>p. (B2f) 



where |X)p, |y)p, |Z)p are the /7-type orbital parts of the Bloch 
amplitudes with wave vector p, which have the same rotation 
and inversion transformation properties as the coordinate sys- 
tem Xp, yp, Zp, defined with respect to the momentum direction 
(i.e., Zp = p/p). The mixing of the orbital wavefunctions and 
the electron spin states in the total angular momentum eigen 
states includes the spin-orbit coupling automatically. This 
spin-orbit coupling is an intrinsic relativistic eff'ect and does 
not reply on whether or not the material has inversion symme- 
try. 

With the convention chosen in Eq. ( |A1[ ), the transformation 
of the Bloch states and spin states are as follows 

|X)p =\X) cos ^p cos (pp-\-\Y) cos ^p sin 0p - |Z) sin 6^, (B3a) 
|7)p = - IX) sin 0p + IF) cos 0p, (B3b) 
|Z)p =\X} sin ^p cos 0p -h |7) sin ^p sin 0p -h |Z) cos ^p, (B3c) 



|T)p = + cos \^) + sin -^e^ 



■ sin \^) + cos ^^^^'^p/^ \i) . 



(B3d) 
(B3e) 



|X), 1 7), |Z) are the orbital Bloch functions which transform 
as X, Y, Z, and | t / |) are the spin Bloch function as the 
eigenstates of cr • Z with eigenvalue ± 1 . With this convention, 
the Berry curvature term has a very simple form as 



-/<7,m1 Vp |j,m)p =i—^Sm'±i,m ^/U ± m)U + m-\-l) 
y2p 



cos ^p yp 

I— . 

sm^p p 



(B4) 



Appendix C: Faraday rotation angle 

For a light with frequency ojq, the polarization density is 

F = €oYj^crXcr,cr'F^^. (CI) 



Then the energy density in the material is 

PE = ^<D(r,0-F(r,0) + ^<B(r,0-H(r,0) 
= (6o6,F -h P) • F* -h c.c. 



(C2) 



where 6^ is the background dielectric constant. Thus the linear 
optical susceptibility is related to the eff'ective Hamiltonian 
through 



cr,o-' 



cr,cr' 



Thus we have 



Xcr,(T' ~^Xcr',cr 



1 ^'Hff 

60 OF^^dF^. • 



(C4) 



The index change due to two circular polarization is respec- 
tively 

Sn± = ^/6r +x±± - ~ -\^--^ ^ = ^\n~^X++^ (C5) 

where n is the material refractive index. The phase delay 
within a propagation length / is then 

d(p+ = ajqC~^lSn+ = 27iA~^lSn+, (C6) 

where A is the light wavelength in vacuum. For a light with 
linear polarization 

x = (-n+-Mi_)/V2, (C7) 
after propagation of the length /, the polarization becomes 

[-n+e-'^"^- -h n-e-'^"^-) / V2 = cos Scp+x -h sin S(P+y. (C8) 
So the Faraday rotation angle is 

(C9) 



Of = S(p+ = :^X++- 



For a light with incident zenith and azimuth angles /3 and y, 
the angles inside the sample JS' and Y are determined by 

nsin/S' = sin/S, (ClOa) 
/ = r, (CI Ob) 

the propagation length through a sample of thickness L is 

l = L/cosfi\ (Cll) 

For a pure spin current and an off'-resonant probe, the suscep- 
tibility is 

/!l = -/y. = i^iqJz + ^2qz • J • z) , (C12) 

46o 

Thus the Faraday rotation for a spin current polarized normal 
to the surface (as in Awschalom's experiment^) is 

^f=^^^ = 0/ 2/cosyg^ sinyg^ cosy 



SeohA cos 13' 



cr,cr' 



2ne{)A^ 

where q = Inn I A has been used. 



sinyScosy, 



(C13) 
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Appendix D: Anisotropic valence band effect 

The anisotropic Luttinger-Kohn matrix of //^^ is 



(E3 P Q ^ 

P* El Q 

e* E.i -P 

2* -P* £_3j 



(Dl) 



where 



£3 = £-3 = ^ [(n + 72)^' - 3r2^z'] , 
£-1 = ^f(ri -r2)^' + 3r2^,2j^ 



£1 



V3y3 
mo 



^zikx iky). 



(D2a) 
(D2b) 

(D2c) 



e = ^[-V3r2(/:^_ ^2)^-2^3^3^^^^]^ (D2d) 
The eigenstates can be in general written as 



(D3) 



By making the transformation U^HUU^a = U^a with 



U 



■t _ 



ro -1 
0010 
0-100 
1000 



(D4) 



we can see U^HU = H*. Thus we get the relation 




FIG. 7: (color online) The geometry of the coordinates. 



U'^a = a. (D5) 



Without loss of generality, we take the coordinate relation 
between the electron and the crystal as 



(100) = sin 0Xp - cos 0yp 

(010) = cos 6 cos 0Xp + cos 6 sin 0yp - sin Oz^ 

(001) = sin 6 cos 0Xp + sin 6 sin 0yp + cos 6zp (D6) 
where (001), (010) and (100) are directions of the three crystal 
axis, and 6 and are the relative direction angles between 
Xp, yp, Zp and the (100), (010), (001) axes. 



The explicit form of Aj^ is 



AX - 



iy - 



smd) -\- 1 



cos + I 



sme/ + 



costycos0 + 



cos ty sm + 



3 («,"')* si 



sm u cos ( 



sin cos 0, 



cos 6, 



(D7a) 
(D7b) 
(D7c) 



and with Eq. (D5 (, they satisfy the relation 



(D8) 
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Appendix E: Second-order nonlinear susceptibility 



With a spin current of the form J = /^XZ + JzZ/L, where Jx is the transverse amplitude, the second-order nonlinear optical 
susceptibility induced by the spin current is 



X 



.(2) 



{oj, , 0^2; a;i + C2) =yx[xXY (-2^3 -2^'^+^',+^',) + TLX (-4^3 - + ^3 - ^5) + YXX (-^4 - ^5 - - ^'^ 
+XYX (-2^3 + ^4 + ^5 - 2^;) + ZYZ (-^1 + ^3 - ^5 - 4^;) + YZZ (^1 - ^3 + ^5 + - + 
+YYY(-4^3-4^;)] (Ela) 
+Jz\ (XYZ - YXZ) (^1 + ^2 + 2^3 + ^4 + 3^5 - ^2 - - 

+ (Yzx - xzY) (6 + 3^3 + ^5 - ri - ^2 - - ^; - 3^;) 

+ (ZXY - ZYX) + 5^3 + ^5 - ^2 - 56 - ^5) 1' (Elb) 



where is derived from by exchanging wi and W2, and 



^1 = 



6 = 



^6, + 2 



^^6, + 2 



\+melmi l-\-me/mi l+m^/m^ l+m^/m^ 



(AO^A^ 
A^ - A^ A 



(AO^^^i (AO^A^ 



i (a^-a^)(a^-a^) 

5 4£^A^A^A^A^ ' 

A^-A^ (a^-A^)(a^ + A^) 



IEfA^A^coi AEpA^A^^A^A^ 



l+nielmh l+nielmh l+nielmi l+nielmi A^ - A^ 



(E2a) 
(E2b) 

(E2c) 

(E2d) 

A^-A^ (a^-A^)(a^ + A^) 



{A^foji 



(A^)2A 



{A^Yoji 



(AO^Ai, EpA^A^coi lEpA^A^A^ 4£^A^A^A^A^ 



(E2e) 



where the factor containing the material dielectric constant 6^ takes into account the difference between the macroscopic external 
field and the microscopic local fieldP 
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